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Abstract 

The excitations of the vortex in Abelian Higgs model with small ratio 
of vector and Higgs particle masses are considered. Three main modes 
encountered in numerical computations are described in detail. They are 
also compared to analytic results obtained recently by Arodz and Hadasz 



> 

00 

O '■ 1 Introduction 

O 

The vortex type solutions of nonlinear field equations are found to be im- 
O^ . portant in many areas of physics. They are useful in considerations concerning 

- ' different phenomena in field theory, cosmology and condensed matter physics 

[gj. However both static and dynamic vortex configurations are known only 
approximately since they are usually described by complicated non-linear differ- 
t^i ential equations. Therefore the numerical methods must be applied to examine 

time-dependent vortex solutions. 

The authors of paper ^ have investigated an excitation of the vortex in the 

Abelian-Higgs model with the mass of the Higgs field much bigger than the mass 

of the vector particle. They considered time dependent, axially symmetric fields 

^ [ with time and radial components of the gauge potential equal to zero. In the 

central part of the vortex the fields were approximated by polynomials. These 
polynomials were continuously matched to the asymptotics describing outer part 
of the vortex. The authors began with finding an approximation to the static 
vortex solution. Next they assumed that the evolution of the Higgs field and the 
azimuthal component of the gauge potential is frozen and considered harmonic 
oscillations of the gauge potential component parallel to the vortex axis. Finally 
they found corrections to the Higgs field and the azimuthal component of the 
gauge potential. 

In the present paper we would like to investigate described above excitation 
of the vortex with the help of numerical computations. We have computed 
evolution of axially symmetric fields with appropriatelly chosen initial data. The 

•Institute of Physics, Jagellonian University, 30-064 Krakow, Reymonta 4, Poland 
^Institute of Physics, Pedagogical University, Podchor§izych 2, 30-084 Krakow, Poland 



approximation used in ||j was simple, and it led to the prediction of oscillations 
of the Higgs field and of the azimuthal component of the gauge potential with the 
frequency equal to the doubled frequency of the oscilations of the gauge potential 
component parallel to the vortex axis. The solutions we have obtained are more 
complicated. Apart from the predicted mode, they contain also oscillations 
corresponding to other vortex modes. 

Our paper is organized as follows. In Sec. 2 we introduce an axially sym- 
metric ansatz and transform the equations to the form which is more convenient 
in numerical computations. Sec. 3 has rather technical character: containes de- 
scription of our computations. The results are presented and discussed in Sec. 
4. Finally Sec. 5 involves some general remarks and conclusions summarizing 
our paper. 

2 Abelian Higgs Model 

The Abelian Higgs model is described by the following Euler-Lagrange equations 

(a, + iqA^w + ^qA")-^ + \h\ * ? -^) = 0, (1) 

^^F^"' = iq ($*9''$ - $9''$*) - 2q^A'' | $ 1^ . (2) 

Our notation follows the paper |1|. Here $ is a complex scalar field, A^, 
is U {!) gauge field, m, q, and A are positive constants. The signature of the 
metric tensor is (-f , — , — , — ) . 

We restrict our considerations to the axially symmetric field configurations 
described by the following Ansatz 



^ = \J^e^'F{t,r), (3) 

^0 = 0, A3 = A{t,r), (4) 

A, = ^sme{l~x(t,r)), (5) 

qr 

A2^~- cos0 1-xt,r , 6 

qr 

where r = ^J2m? {{x^Y + {x'^Y), 9 = arctan (x^/x^), t = \/2rna;°. The fields 
have to be non - singular on the x^ axis. This requirement implies 

F(t,0)) = 0, x(i,0) = l, |-F(r = 0,t) = 0. (7) 

The axially symmetric ansatz @-@ applied to the equations (Q)-(||) simplifies 
them to the form 



F = F" + If' - f^x' + ^aA + i (F - F3) , (8) 

X^Xii- -XI - li^F^X, (9) 

r 

A = A" + -A' - K^F'^A, (10) 

r 

where dot denotes the time derivative and prime the derivative with respect to 
the variable r. The dimensionless parameter k = ^2q^ jX is equal to the ratio 
of the vector particle mass and the mass of the Higgs particle. In this paper 
we restrict ourselves to small values of k (0.05, 0.1, 0.2). We also assume that 
qjm = 1. This does not Hmit the generality of our considerations since it can 
be achieved by rescahng the field A. 

Let us note that all field configurations of the form (Q)-@ with finite energy 
per unit of length in the x^ direction have unit topological charge. If in addition 
we assume that the functions / and x ^^re time independent solutions of the 
equations @-(||) with A3 = A = then we obtain the well known Nielsen- 
Olesen vortex. This static solution was examined in many papers. However its 
exact analytical form is not known. Only approximate methods (both analytical 
and numerical) have been worked out and can be adopted as the starting point to 
the precise numerical calculations. These are neccessary as the time dependent 
analysis is very subtle. 

The equations (p|)-(|l0|) are singular in r = 0. Therefore we apply to them 
the following transformation which removes the singularity 

r2 
^=^^> (11) 



/(i,a;) = ^^^M, (12) 

r 

Mi,a;) = ^M^, (13) 

a(i,x) = A(i,r). (14) 

Thus we obtain the equations of the form 




(2/,, + r2;i2^a2)/ + i/(l-r2/2), (15) 



h = (^-^^^] h" + {—^ + ^i^l h' - n-'f (1 + r^h) , (16) 




a' - K^r^fa. (17) 

Here prime denotes the derivative with respect to x, r = (x+ \Jx (4 + x) j /2. 

Some disadvantage of the transformation (O) is that the neighbourhood of the 
point r = is mapped into a very small x interval (dx/dr |r=o= 0). Therefore 
the detailed solution in this region has to be examined separately. 

3 The method 

In this paper we investigate the oscillations of the vortex caused by the ex- 
citations of the field A. The first step is to have the precise solutions for the 
functions / and % in the static vortex configuration. In ^ these fields in the cen- 
tral part of the vortex were approximated by polynomials. These polynomials 
were continuously matched to the appropriate asymptotics describing external 
part of the vortex. There also exist more precise numerical approximations of 
the static vortex |^ but their accuracy is not sufficient too. Therefore we have 
started with one of these approximate methods and then in order to improve 
the solutions we have applied the average method, that is we have examined 
the time evolution of the vortex alone averaging the functions / and x around 
their equilibrium state. 

The next step is to excite the vortex. We have done this by adding a non- 
zero field A. Following ||l| we have used a configuration which gives harmonic 
oscillations of this field while the vortex fields are frozen. Substituting A{t, r) = 
a{r) cos cot into (ttOh one obtains 

a" + -a' - K^F^a -I- w^a = 0. (18) 

r 

As was noticed in [Q this equation has the form of the one-dimensional 
Schrodinger equation and possesses at least one bound state. Rewriting it in 
the variable x we obtain 




-(KV/2-a;2)a = 0. (19) 

On the uniform grid with the spacing A eq. ( [l9| ) leads to the following finite 
difference equations 

— {di-aa)+Lj^do^O, (20) 






[l + rY \A 



,2 



- {n^^f - u^) flfe = 0, fc = 1, 2, ... (21) 

For small values of the parameter k the frequency uj is approximately equal to 
K ly . Therefore we can solve the above difference equations (|2(]| ) , ( pi|) where w 
is replaced by k. We choose arbitrary clq and succesively determine oi, 02, ... . 
The graphes of the function A{r) — a(r^/(l + r)) obtained from formulae (po|), 
( pi| ) with A = 0.01 and the approximation of the function A obtained in the 
paper [Q for k = 0.1 are ploted in Fig 1. 

Now we are going to investigate the evolution of the excited vortex using 
numerical computations. We solve the following finite difference equations 

f]+i,k = 2/j^fc -/j_i,fe + f - j [t^—^] {f].k+i-2fj,k + fj.k-i) 

i)'A(2|^ + ^) (/,,,,-/,,) (22) 



2 



iV I (2 + r)r 



2 



'ij+i,fc — 2/ij,fe — hj^i^k + { - ] I 7:7- — "2" I [hjM+i — '^■hj,k + hj^k-i) 



V \(i 



r 



^\\ J^ + ^f 



4 J ^i^^TTo^ + ^TTo^)^^^'--'^^-^'^) 



(^)'A^-^(i+r^%.)/i:.., 



(23) 



2 



2 



,3V / (2 + r)r , , ^ , 

aj+i,fc = 2aj.k - aj-iM + 1 4 ) I J Z2 I (Oj.fc+i - 2aj.k + aj,k~i) 

^\\ /^o(2 + 0' 4 + r ^ 



4 J ^i'aT^y^^aTo'j^"^^'^"^""^'''^ 

ly A'.\\,kflk, (24) 

where A is the space between the grid points in the x-direction; the space 
between time slices is 3 A/4. As the initial values for the fields / and h at t — 



we take found earlier functions corresponding to the static vortex solution. The 
initial values of the field a are computed from the formulae (pO|), ( pi| ) with 
ao = 0.1. We also assume that initially the derivatives with respect to time of 
all fields are equal to zero. In order to remove disturbances of the solution caused 
by the right boundary of our grid we have performed computations on the larger 
area and then we removed the part of the grid which could be disturbed. 

4 Numerical results 

We have computed the evolution of the vortex for < i < 600 and k, = 
0.05,0.1,0.2 taking A = 0.01. The field A oscillates with the frequency a; w k. 
Fig.2-Fig.7 present the evolution of the functions 5fx{t) = f{t,x) — f{0,x), 
Shx{t) = h{t,x) — h{0,x) (i.e. the differences of the time-dependent and static 
values of the fields / and h respectively) at the points x = 5 and x = 160. 
We have chosen these values of the variable x to compare the oscillations near 
the vortex core with those far from it. The authors of the paper ||l|| predict 
oscillations of the field / and x with the frequency ui — 2k but they have only 
analysed the modes forced by the oscillations of the field A. The numerical 
results we have obtained are more complicated and involve also other vortex 
modes. These modes will be presented below in more detail. It is also interest- 
ing that while the oscillations of the Higgs field / start at once for all values of 
X the oscillations of the field h are delayed for points far from the vortex centre. 
They begin at the vortex core and then propagate outside. This fact is also 
illustrated in Fig.8-Fig.l0 which show the plots of the functions 6hx for fixed 
time values t = 100, 200, 300. One can see that in the region r > t the function 
6hx is almost equal to zero. It starts to oscillate only when the disturbance of 
the field h arrives from the vortex centre. 

Let us now present the frequencies of the vortex oscillations in the neighbour- 
hood of its equihbrium state. We have computed the fourier transformations 
for the functions 6fx and Shx for three different values of k mentioned above 
eg. K = 0.05, 0.1, 0.2. These calculations have been performed for two distances 
from the vortex core. As before we have choosen points x = 5 and 160. We 
have presented our results for the Fourier transformations in Fig.ll-Fig.l6. The 
most important conclusions can be summarized as follows. The main frequency 
of the vortex oscillation is equal to 2k. The appropriate peak can be easily 
observed in each of the figures Fig.ll-Fig.16. These oscillations are explained 
in detail in ||l| on the basis of the approximate analytic calculations. They are 
forced by the corresponding oscillations of the field A with frequency k. But 
the other vortex modes are also excited. Let us stress that their frequencies 
do not depend on the parameter k. The approximate numerical values of these 
frequencies are correspondingly equal to the mass of the Higgs particle (1.0 in 
our case) and 90% of this mass (that is 0.9). Let us note that these values seem 
to be independent from the field a and are two points in the whole spectrum of 
frequencies of the vortex oscillations computed in |Q| . However it is interesting 
that the lower frequency was not detected for the field / far from the vortex 



core (Fig. 12). 

5 Ending Remarks 

The excitation of the vortex we have investigated was obtained by choosing 
some particular initial data for the component A^ of the gauge field parallel to 
the vortex axis. In the static vortex solution this component of the gauge field is 
equal to zero; in our case it oscillates with the frequency approximately equal to 
k: the ratio of the vector boson and Higgs masses. The other fields also oscillate. 
Performing the Fourier transformation we found that the frequencies of these 
oscillations are grouped near 2k, 0.9, 1.0. The first value was predicted in ||l|| and 
is caused by the term proportional to A^ in (g). The oscillations with the second 
frequency appear mainly on the vortex core and are probably concerned with 
the bound state of the static vortex solution |Q]. The third frequency simply 
equals to the mass of the Higgs particle. However it should be noted that our 
computations include about ten periods of the field A oscillations and therefore 
the question concerning the long time behaviour of the system is still open. The 
field configuration we have considered was constructed in such a way that the 
field A starts its oscillations at once in the whole space. These oscillations excite 
the Higgs field for all values of the spatial coordinate. However the motion of 
the azimuthal component of the gauge field looks differently. It begins at the 
vortex centre and then propagates outside. 
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